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ABSTRACT 

We construct the equation of state (EoS) for neutron stars explicitly including 
hyperons and quarks. Using the quark-meson coupling model with relativistic 
Hartree-Fock approximation, the EoS for hadronic matter is derived by taking 
into account the strange (a* and <f>) mesons as well as the light non-strange (cr, 
oj, 7r and p) mesons. Relevant coupling constants are determined to reproduce 
the experimental data of nuclear matter and hypernuclei in SU(3) flavor symme¬ 
try. For quark matter, we employ the MIT bag model with one-gluon-exchange 
interaction, and Gibbs criteria for chemical equilibrium in the phase transition 
from hadrons to quarks. We find that the strange vector ( 0 ) meson and the Fock 
contribution make the hadronic EoS stiff, and that the maximum mass of a neu¬ 
tron star can be consistent with the observed mass of heavy neutron stars even 
if the coexistence of hadrons and quarks takes place in the core. However, in the 
present calculation the transition to pure quark matter does not occur in stable 
neutron stars. Furthermore, the lower bound of the critical chemical potential of 
the quark-hadron transition at zero temperature turns out to be around 1.5 GeV 
in order to be consistent with the recent observed neutron star data. 

Subject headings: dense matter — elementary particles — equation of state - 
stars: neutron 


2 


Introduction 


White dwarfs, neutron stars, and black holes, which are collectively referred to as 
compact stars, are consequent remnants of core-collapsing supernovae explosions. Neu¬ 
tron stars may especially be believed to be cosmological laboratories for nuclear matter 
at extremely low temperature and high density, because the central density of neutron 
stars can reach several times higher than the normal nuclear density. Thus, pulsar ob¬ 
servatio ns can provide some constraints on the equation of state (EoS) for dense nuclear 


matter fl Weber! 11999c iGlendennind 120021; lLattimer fe Prakashl 120061 1 . Of particular interest 


is the possib ility of exotic degrees of freedom in the core of a neutron star, such as hy- 
peron s (Ys) f Glen denning Rr. Moszkowskil 1991 : Schaffner fe Mishustin 1996 : Nishizaki et al. 


of boson-like matter (Takatsuka et al 


199811 and/or dark matter (jPerez-Garcia et all 2010) 


2002 1. quark matter ( Itoh 197ol: Witten! 1984 : Weberl 2004 h_ some unusual condensations 


1978l~~ Tatsnmi 19881: Idendenning fe Schaffner-Bielich 


Thanks to recent advances in astrophysical observations, we can obtain some precise 
information on the properties of neutron stars. In particular, the discovery of massive neutron 
stars, PSR J16 14-2230 with 1.97±0.04 Mm ( Demorest et al. 2010 1 and PSR J0348+0432 with 
2.O1±O.O4M 0 (iAntoniadis et all 20131 ). sets a strong constraint on the EoS for dense matter. 
Meanwhile, at present, there have been many theoretical studies of the nuclear EoS based 
on many-body theories. However, it is quite difficult to explain the heavy neutron stars by 
the EoS which have been calculated so far, if hyperons are supposed to exist in the core of 
neutron star, because the degrees of freedom of hyperons make the EoS very soft, and thus, 
the possible maximum mass of a neutron star is considerably reduced. 

In order to solve this discrepancy between the observations and theories, which is the 
so-called hyperon puzzle, several useful approaches have been proposed in the last few years. 


model ( 

Mivatsu et al. 

2012; 

Katavama et al. 

2012; 

Whitten 

)urv et al. 

2012) and the Dirac- 

Brueckner-Hartree-Fock approach ( 

Katavama & Saito 

2013, 

2015 

) show remarkable results 


to settle the hyperon puzzle. Their calculations indicate that the EoS with hyperons keeps 
stiffness even at high densities due to the suppression of hyperon production and a up¬ 
ward shift of the density at hyperon appearance. In the relativistic mean-field (RMF) [or 
relativistic Hartree (RH)] calculations, the extension of SU(6) spin-flavor symmetry based 
on the quark model to SU(3) flavor symmetry with strange mesons is also helpful to clar¬ 
ify the problem in determining the couplings of the mesons to the octet baryons, because 
the strange vector (<t>) m e son plays an importa n t role in supporting ma s sive neutron stars 


flWeissenborn et al. 


2011 


Mivatsu et al.l l2013al: iLopes fe Menezesl I2014J: I Jiang et al.l 12012 


Co lucc i fe Sedra kian 20131). It is also worth studying to understand the effect of hyperons in 
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dense matter, introducing tl 

Le interaction terms among various mesons ( 

Bednarek et al. 

2012; 

Sulaksono & Agrawal 

2012; 

Tsubakihara & Ohnishi 2013h or using quantum Monte Carlo 


calculations flLonardoni et al.ll2015h . Furthermore, t he modified f (R) grav i ty may be another 


cand idate for the solution of the hyperon puzzle (jCheoun et al.l 12013d: lAstashenok et al. 
2014h . 


Although the phase transition from hadrons to quarks at high temperature and low 
baryon chemical potential is known to be crossover as explored by the first-principle lat¬ 
tice QCD simulation flAoki et al.l 120061 ). not only the order of the phase transition at zero 
temperature but also the existence of a critical end point in the QCD phase diagram are 
still unknown. In dense matter such as a neutron star, the possibility of quark matter as 
well as hyperons may be expected to appear in the core and to influence a number of in¬ 
teresting astrophysical phenomena. In general, the EoS for neutron stars with hadrons, 
leptons, and quar ks, which is so-called hyb r id stars, is conside red by assuming a first-order 
phase transition (jClendennind Il992l. l200ll; lAlford et al.l 120131 ) . There have recently been 


many papers where the maximum mass of neutron stars with deconhn ed q uarks ca n e xcee d 


2M m if the vector interaction between quarks is strong enough ( 

Bonanno & Sedrakian 

2012; 

Lenzi & Lugones f 

>012|; I 

jogoteta et al. 2 

Old; 

Logoteta & Bombaci 

2013; 

Orsaria et al. 

2013, 

2014; 

Klahn et al. 

2013|; 

Yasutake et al 

2014). In addition, a 

smooth crossover based on 


the percolation picture to obtain the quark-hadron transition at zero temperature has been 
studied, and the E oS can also sustain th e neutron-star mass of 2 M e using strongly interact¬ 
ing quark matter (Masuda et al. 2013allbl h As explained in Refs. ( Lenzi &; Lugones 2012 


Orsaria et al. 2014), without the strong vector interaction between quarks, it would be hard 
to explain the observation of compact stars with mass greater than around 2 M & . 

In this paper, we construct the EoS for neutro n stars w hich ca n sat isfy the 2AR. con¬ 


strain t from the recent astrophysical observations (IDemorest et al.l l2010l; lAntoniadis et al. 


201 31). even when the exotic possibility such as hyperons and quarks are considered in the 
core. Then, we study the properties of neutron stars, for instance, the mass, radius, and 
particle fractions, especially by focusing our mind on the effect of strangeness in neutron 
stars. 

For uniform hadronic matter, we employ our previous EoS which have been calcu- 
l ated us ing the chiral quark-meson coupling (CQMC) model within RHF approximation 
fiMivatsu et al.l 20 1 3a ). The CQMC model is an extended version of the QMC model, in 
which the quark-quark hyperhne interactions caused by the one-gluon and pion exchanges 
are incl uded. Such hyperhne interactions play an impotent role in the baryon spectra in 
matter ( Nagai et al. 2008 : Mivatsu fe Saito 2009 : Saito 2010 b 


In the QMC model, the quark mass in nuclear matter is reduced from the value in 
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vacuum because of the condensed scalar (a and a*) fields depending on the nuclear density. 
The decrease of the quark mass then leads to the variation of baryon internal structur es at 
the quark level. Such effects are considered self-consistently in the QMC model (IGuichon 
19881; Saito & Thomasl 1994aj). In fact, the evidence for the medium modification of nucleon 
(IV) structure in a nucleus has been observed in polarization transfer measurements in the 
quasi-elastic (e, e'p) reaction at the Thomas Jefferson Laboratory, and the result supports 
the prediction of the QMC model (IBrooks et al.l 1201 lh . This model has been successfully 
applied to the European Muon Collaboration effect (IGeesaman et al.lll995l: ISaito fe Thomas 
1994bl: ICloet et al.l 120061) and to exp l ain various properties of finit e nuclei as well as infinite 


nuclear matter (IGuichon et al.lIl996l: ISaito et al. 


19961.119971.120071) . This approach can a lso 


be useful to evaluate the density dependence of various form factors (jCheoun et al. 2013a.Ij). 


In addition to the CQMC model, relevant coupling constants are determined so as to 
reproduce the experimental data of nuclear matter a nd hypernuclei in SU(3) fl avor symmetry, 
includ i ng the hidden strange (a* and 0 ) mesons (IWeissenborn et al.l 1201 ll: iMivatsu et al. 
2013a : Lopes fe Menezes 2014 ). 


In order to study how quark matter affects the neutron-star properties, we adopt 


trary to the recent calculations with it (Bonanno & Sedrakian 

2012- 

p — ? p 

Lenzi & Lugones 

2012; 

Logoteta & Bombaci 

2013: Logoteta et al. 2013: Orsaria et al. 

2013, 

2014; Klahn et al. 

2013; 

Yasutake et al. 2014) 

. Furthermore, we assume the first-order phase transition from hadrons 


to quarks under Gibbs criteria flGlendenuiugill992l. 120011) . 


This paper is organized as follows. In Section [2l a brief review ofjdny formalism for 
Hartree-Fock calculation based on quantum hadrodynamics (QHD) flSerot fc Walecka 198 4) 
is presented. The description of quark matter with one-gluon exchange effect using the MIT 
bag model and the phase transition from hadronic matter to quark matter are explained in 
Section [3] and [4j respectively. Numerical results and discussions are addressed in Section [5] 
Finally, we present a summary in Section [6] 


2. Description of hadronic matter 


We present the formulations for des cribing uniform hadronic matter. In Quantum 
Hadrodynamics (QHD) flSerot fe Waleckal 1 9841 ). the baryons are treated as point-like ob¬ 
jects, and interact via the exchanges of mesons. On the other hand, we want to include the 
effect o f baryon-structure variation in matter using the chiral quark-meson coupling (CQMC) 
model fjNagai et al.ll2008l: IMivatsu &; Saitoll2009l: ISaitoll2010l ). In this paper, the Lagrangian 



































































































































density for hadronic matter is thus chosen to be 



( 1 ) 


where 


^ b (*7a»^ ~ M b ) iI>b, 


( 2 ) 


B 


with i^b being the baryon field and Mb being the baryon mass in a vacuum. The sum B 
runs over the octet baryons: proton (p), neutron (n), A, £ +0_ , and S°~. For the free baryon 
masses, we take Mjy = 939 MeV, M\ = 1116 MeV, Ms = 1193 MeV, and M= = 1318 MeV, 
respectively. Lepton Lagrangian is introduced in Section [H 

In the present calculation, we study the effects of direct and exchange contributions on 
hadronic matter through not only the exchanges of non-strange mesons (cr, cu, 7r, and p) but 
also those of strange mesons (cr* and 0). Thus, the meson term reads 



( 3 ) 


with 


( 4 ) 

( 5 ) 

( 6 ) 



where the meson masses are respectively chosen as m a = 550 MeV, m a * = 975 MeV, 
m u = 783 MeV, rri^ = 1020 MeV, m p = 770 MeV, and m n = 138 MeV. 

The interaction Lagrangian is given by 



where the common mass scale, A4, is taken to be the free nucleon mass, and Ib is the isospin 
matrix for baryon B. The cr-, a*-, u 0-, p-, tt-B coupling constants are respectively denoted 








6 


by g<rB(cr), g<r*B(cr*), g^B, g^B, g P B and f wB , while f uB , Ub and f pB are the tensor coupling 
constants for the vector mesons. In the CQMC model, the coupling constants, g^sic r ) and 
g**B(cr*), have the scalar-field dependence which reflects the variation in the internal (quark) 
structure of baryons in matter, and they are caused by the attractive interactions due to 
the a and a* exchanges. For simplicity, we adopt the following simple para metrizations for 
those coupling constants flMivatsu et al.ll2013ai I20l4 iTsushima et al.l 1 199811 : 


g<jB{<y ) — S'crs&s 


9a*B (&*) — 9a*B^>l 


B 


1 - Y (fl 1<tNO) 

1 - Y (<wO 


( 8 ) 

(9) 


where g a n and g a *a are respectively the a-N and cr*-A coupling constants at zero density. 
The effect of the variation of baryon structure at the quark level can be described with the 
parameters a# and a' B . In addition, the extra parameters , bn and are necessary to express 


the effect of hyp erfine interaction between two quarks fjNagai et al.l 12008c iMivatsu fc Saito 
20091 : Saito 2010). The couplings in the CQMC model are invariant under Lorentz transfor¬ 
mation because they are functions of the scalar fields, and the values of the four parameters 
given in Equations (J8]) and flSl) are tabulated in Table [U If we set a# = 0 and b B = 1, g a B(&) 
becomes identical to the a-B coupling constant in QHD. This is also true of the coupling 

9a*B(<7*). 

In mean-field approximation, the meson fields are replaced by the constant mean-field 
values: a, a*, u ;, 0 , and p (the p° field). The mean-field value of the pion vanishes in 
the relativistic Hartree (RH) calculation, while the pion effect should be included in the 
relativistic Hartree-Fock (RHF) calculation, where the exchange contribution as well as the 
direct one are taken into account. To sum up all orders of the tadpole (Hartree) and exchange 
(Fock) diagrams in the baryon Green’s function, Gb, we use the Dyson’s equation 


a B (k) = G° B (k) + G%(k)S B (k)G B (.k), 


( 10 ) 


where kJ l is the four momentum of baryon, Eg is the baryon self-energy, and G° B is the 
Gree n’s function in free space. The baryon self-energy in matter is generally written as 
( Serot fe Walecka 1984 1 


s B (k) = S’ B (k) - 7„E» (fc) + (7 ■ k)S’ B (k), 


( 11 ) 


with k being the unit vector along the (three) momentum k and being the scalar 

part (the time component of the vector part) [the space component of the vector part] of 
the self-energy. Therefore, the effective baryon mass, momentum, and energy in matter are 
respectively defined by including the self-energy in matter as follows flBouvssv et ah 119871: 
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Mivatsu et al. 

2012; 

Katavama et al. 

2012) 


M‘ B (k) = M b + E %(k), 
kf = (fcg>, k%) = (fc° + Eg(fc), k + k^ B (k)), 
E B (k) = [kf + M# (*)] 


1/2 


( 12 ) 

(13) 

(14) 


The baryon sel f -energies in Equation (TTlT) are then calculated by (jBouvssv et al.lll987l; 
Mivatsu et al. 2012 ; Katavama et al. 2012 1 


= -9crB(cr)o ~ g<r*B(v*)v* 


E 


i^BB ') 


2 rkf 


(47 l) 2 k Jn 

B',i y ’ JU 


dq q 


M B'(q) p n \ . g*B' D , k 4 

[ Eg,(q) Bl ^ q) + 2E* B ,(q) 


— —g^B^ — sv>,B0 — g P B{iB)zP — ^2, ^ BB ^ ^ 


B'i 


( 47r ) 2fc A 


dq qAi(k,q), 


£b(*0 = E 


(Ibb 1 ) 


2 /-fc_F 


B r; (47r)' 2 fc 7 0 


dg g 




.Eg/iq) 


' / t \ Ad d/ (g) / 7 \ 

T~sCi(k, g) + , iA(fc, d) 




(15) 

(16) 
(17) 


where kp B is the Fermi momentum for baryon B, and the factor, I BB i, is the isospin weight 
at the meson-Bd?' vertex in the Fock diagram. 

In order to include the effect of the finite size of baryons, a form factor at each interaction 
vertex shoul d be introdu ced, and, in the present calculation, we employ a dipole-type form 
factor (Katavama et al. 2012): 


Hp 2 ) = 


(1-P 2 /AD 


2 i 


(18) 


where is the (four) momentum transfer, A* is a cutoff parameter, and i specifies the 
interaction vertex (see the 1st column of Table [2]) . In the interaction Lagrangian density, 
we then replace all coupling constants with those multiplied by the form factor. In fact, the 
effect of form factor can not be seen at the Hartree level, because the momentum transfer 
in the meson exchange between two baryons vanishes. However, because the exchanged 
momentum can be finite in the Fock term, it may become significant as the density increases. 
In addition, the retardation effect in the Fock ter m are ignored, since it gives at most a 


few percent contribution to the baryon self-energy (ISerot fe Waleckal 1 19841 : iKatavama et al. 


20121; IWhittenburv et al.l l2Q12h . The functions A if B i5 Cj, and Di in Eq uations f[T5l) - ffTTl) 


are explicitly given in Table [2j in which the following functions are used (iKatavama et al. 

















































2012 ): 


Qi(k,q) 


®i(k,q) 


vi{k,q ) 

n i(k,q) 

r i(k,q) 




' M?(k,q) L7(k,qy 
Mr(k,q) Lf(k, q) _ 


( m i - A2 i) n K(k,q) 


n= 1 


4 ^ [( fc2 + 9 2 + m 2 ) Oi(k, q) - A 8 N?(k, q)] , 

( k 2 + g 2 - m- / 2 ) $*(/;:, q) - kqQ t (k, q) + fg), 
[k 2 + g 2 ) $*(£;, g) - kq<di(k, q) + g), 


[fc0i(A:,g) - 2g$i(fc,g)], 


where 

A 8 

^i{k, q) = ( k , g) + ( k 2 + g 2 + A 2 ) Nf(k, g)] , 

Lf(k,q) = A 2 + (fc±g) 2 , 

Mf(k,q) = m 2 + (fc ± g) 2 , 

JVT(*.4) = ([L+(M)P - [£,-(*. 4)]"”) . 

We note that the effect of form factor vanishes in the limit A* —> oo. 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


(24) 

(25) 

(26) 

(27) 


By solving the Euler-Lagrange equations for the meson fields, as usually adopted in the 
relativistic mean-field approach, the me 
and (USD are given by 


a = 


B 

E 

B 


CO = 


B 


B 


?=£ 


•field values of cr, cr*, 

co, (j), and p in Equations (1T5|) 

^C B (a)p s B , 

(28) 

rn a* 

(29) 

- 2 "ds, 

ml 

(30) 

g<t>B 

—yds, 

(31) 

—j(Ib)3PB, 

(32) 


B 


where the scalar density, p s B , and the baryon number density, ps, read 


Pb ~ 

Pb = 


27 b + 1 
27r 2 

2.7b + 1 

2 vr 2 




dk k 


E%(k) ’ 


" B dkk 2 = 2Jb 1 h 3 


67T 2 


Fb ’ 


(33) 

(34) 
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with J B being the spin degeneracy factor of baryon B. In Equations 


and 


C B {a) 


and C B (a*) are respectively the scalar polarizabilities at the a-B and a*-B interactions. In 


the CQMC model, they can be exp ressed by the following parametrizations (IMivatsu et al. 
2013a . 2014 : Tsushima et al. 1998h : 


C B (a) = b B [1 - a B (g<rNcr)] , (35) 

C' B ((j*) = b' B [1 — a ' B {da* vd*)], (36) 

where the parameters, a B , b B , a' B , and b' B , take the same values as in Equations (JSD and ([UJ). 


With a self-consistent calculation of the baryon self-energies in Equations (fT5]) - ([TT|) . the 
energy density for hadronic matter, which includes the baryon and meson contributions, can 
be expressed as 





T B (k) + \v B (k ) , 


(37) 


with 


T B (k) 

V B (k) 


M B M* B (k ) + kk* B 

Wk) ’ 

Em 


Then, the pressure for hadronic matter is given by 


Y?m. 


Ph — n B 


d 


tH 


dn B \n B 


where the total baryon number density is defined by n B = p B . 


(38) 

(39) 


(40) 


3. Quark matter description 


We briefly present the desc ription of uniform quark matter. The thermodynamic poten¬ 
tial can be simply expresse d by ( Freedman fc McLerran 1978 : Farhi fc Jaffelll984 : Alcock et al. 
1986 : Haensel et al. 1986 1 

n = ^n q + B, (4i) 


with the quark term, Q q , and the bag constant, B. The sum q runs over three-flavor quarks 
(u, d and s), and the quark thermodynamic potential is given by a sum of the kinetic term 
and the interaction term due to the one-gluon exchange (OGE) ({F reedman & McL erran 
1978l l: 

9q (2 jq + 1) 


kl q = 


24tt 2 


F (fji q , m q ) - — G {p q , m q ) 

7r 


(42) 
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where the color weight for quark species, g q , is equal to 3, J q is the the spin degeneracy 
factor, and a c is the QCD fine structure constant. In addition, the first-order correction due 
to the gluon interaction is included in the second term of Equation (I42[) . The functions F 
and G are written as 


F (p q , m q ) = fj. q yJfj.l -m 2 q (u 
G(n q ,m q ) = 3 


5 2 \ , 3 4l 
-ml + -m ? In 


Pq *h \fp\ 


mi 


m n 


111 


p q + Vp 2 q ~ m ‘i 


m n 


2 (p 2 q -m 2 q )\ 


(43) 

(44) 


with fig being the quark chemical potential and rn q being the current quark mass. We here 
take m u (rf) = 5 MeV and m s = 150 MeV. 


The quark number density is related to fl q via p q = —dQ q /dp q , and the baryon density 
and the charge density in quark matter are given by tib = | Pq an d n c = | Pu~\ (Pd + p s ), 
respectively. The energy density and pressure for quark matter are then written as 


£q — ^ (Dg + flgPg) + B, 

q 

(45) 


(46) 


g 


Instead of hadronic matter, quark matter is expected to exist at very high densities 
such as the center of a neutron star. However, the exact value of the transition density from 
hadron phase to quark one is still unknown. In the present calculation, in order to study 
the effect of quark matter on neutron stars, we introduce a density-dependent b ag c onstant, 
which is assumed to be given by a gaussian parametrization (Burgio et al. 2002a.b) 


B (n B ) = Boo + (-Bo - Boo) exp 


-i>(^ 

\n 0 


(47) 


with no and /? being the normal nuclear density and a parameter, respectively. Since quark 
matter may appear at very high densities and the final results do not sensitively depe nd on 

Mr 


and Bqq, we here fi x B 0 = 400 MeV fm 3 and B^ = 50 MeV fm 3 (IMaieron et al.ll2004 


Maruvama et al. 20071 ). Due to the density dependence, the bag constant in Equation 
becomes B (n B ), and that in Equation (1461) should be replaced as 

dB ( n B ) 


B —>■ n B - 


dn 


B 


B (n s ). 


(48) 


The coupling between a quark and gluon is scale-dependent, and the lowest-order cou¬ 
pling at momentum transfer Q 2 is given by a c (Q 2 ) = 12n/ [(33 — 2Nf) In (<3 2 /Aq CD )] 
























-11 - 


with Nf being the number of quark flavors and Ag co ~ 200 MeV. In practice, it can be 
parametrized in a convenient form (jCapstick & Isgur 1986j) 

a c ( Q 2 ) = 0.25e" Q2 + 0.15e- Q2/1 ° + 0.20e- Q2/10 °°, (49) 

where Q 2 is in GeV 2 . The coupling in Equation (l49li shows a c = 0.6 at zero momentum 
transfer, and the momentum trans fer is replaced by the average of quark chemical potentials, 
/i 2 = | (/i 2 + p 2 d + p 2 ) (Freedman & McLerran 1978). 


4. Neutron-star matter and phase transition 


In order to calculate the properties of neutron-star matter, the charge neutrality and j3 
equilibrium under weak processes are imposed. The leptons must be introduced to achieve 
these conditions, and the Lagrangian density can be written as 

- me) if) t , (50) 

i 


with 'ipe being the lepton held and me being its mass. The sum £ is for e and p . The 
lepton energy density, pressure, and number density are then given by 


Q = 


Pt = 


E 


2J t + 1 


rkF, 


E 


2vr 2 J 0 

2Je + 1 


dk k 2 \ k 2 + 


m 


ii 


K F e 


Pe = 


3 ^ 2vr 2 

t 

2Je + 1 

27T 2 


dk 


k A 


y/k? + mj 


Ft 


dk k 2 = 


2Je + 1 

67T 2 


k F(’ 


(51) 


(52) 


(53) 


where Jp is the spin degeneracy factor of lepton t. The total energy density and pressure 
for hadronic (quark) matter are given by the sum of the hadron (quark) and lepton parts, 
namely e = t p(o\ + £/ and P = Phiq\ + Pe- Furthermore, the condition of /3 equilibrium is 
expressed as f Glendenning 1992 . l200ll : Maruvama et al. 2007 1 


l^n hA k u T 2p, d , (^^) 

pn T Pe PY,— p3~ — Pd — Psi (55) 

Pn Pe = pp = PY+i (56) 

pe pp,: (57) 


with the chemical potentials for baryons, quarks, and leptons. We note that the properties 
of neutron-star matter are generally characterized by two independent chemical potentials 
related to baryon number density and total charge density. 































12 


In order to describe t he coexiste nce of h adron s and quarks, we impose Gibbs criterion 
for chemical equilibrium ( Glendenn i na 1992, 2001). Under Gibbs criterion, in the mixed 


phase (MP), pressure in the hadron phase (HP) must balance with that in the quark phase 
(QP) to ensure mechanical stability as follows, 


-pHP (pm pe) PqP {.pni he) ■ (58) 

In the MP, where the condition, Equation (158]) . is satisfied, the charge neutrality can be 
expressed as 

(1 - X>c P + X n W = 0, (59) 

with n^ p (n^ P ) beign the charge density in the HP (QP), and x being a volume fraction. 
The total energy density and baryon number density in the MP are then given by 


£mp — (1 — x)eHP + yeqp, 

nT = (1 - X)"S P + xnf. 


(60) 

(61) 


5. Numerical results 


5.1. Coupling constants and matter properties 


Firstly, we consider the iso-symmetric nuclear matter around the normal nuclear density, 
no = 0.155 firm 3 , in which region hyperons do not appear. Then, the couplings for g a N, g u n, 
and g^N are determined so as to fit the saturation energy (—16.1 MeV) at n 0 . We note that, 
in SU(6) symmetry, the (j) meson does not couple to the nucleon (g^N = 0), but that it does 
couple to the nucleon in SU(3) symmetry through the following relation: 


9 <t>N — 


\^3z — tan 9 V 


~9u jNi 


(62) 


1 + \/3z tan 9 V 

where the mixing angle and the ratio of a coupling for the octet state to o ne for the singlet 
state are respectiv ely chosen to be 0 V = 37.50° and z = 0.1949 (jMivatsu et al.ll2013al. I2014J; 


Rijken et al.l 20101). The coupling constants, g a n, g^N, g<pN, and g P N, are shown in Table [3] 
For comparison, we also show the results for the case in SU(6) spin-flavor symmetry. The 
p-N c oupling constant, q„jy, is determined so as to fit the the symmetry energy, E sym = 32.5 
MeV (ITsang et all 2012; Lattimer & Liml 2013). For the tensor coupling constants, we use 
th e ra ti os, J^n/q^n, and f P N/g P N, which have been suggested by the ESC08 model 

(jRi jken et alJl2010l n. Furthermore, the n-N coupling constants, f n n, are also chosen to be 


1 The ESC08 model may at present be the most complete meson-exchange model based on SU(3) flavor 
symmetry, which can describe the Y-N and Y-Y as well as N-N interactions. 
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the value given in the same model. In the present calculation, we assume that the a* meson 
does not couple to the nucleon (g a *N = 0). 

In Table 0] we present the several properties of symmetric nuclear matter at n 0 . The 
symmetry energy, E sym , is defined in terms of the 2nd derivative of the total energy with 
respect to the difference between proton and neutron densities. Then, the slope (curva¬ 
ture) parameter of the symmetry energy, L 0 (A' sym ), is evaluated by th e 1st ( 2n d) deriva- 
with respect to the baryon density f Chen et al. 20091: Agrawal et al. 


tive of E Kvm with res 

pect 

Sulaksono & Agrawal 

2012 ) 


2012 


The incompressibility in the SU(3) case is slightly smaller 
than that in the SU(6) case, and it stands closer to the range of Kq = 240 ± 20 MeV, 
which is derived from the isoscalar giant monopole resonance fjChen et all 20091). The slope 
parameters of the symm etry energy in both case s show the reasona ble values predicted by 
theoretical calculations (IDanielewicz fe Leell2009l; iTsang et al.l 120121) . Moreover, in both of 
the SU(6) and SU(3) symmetries, the 2nd derivative values of the isobaric incompressibility 
coefficient lie well within the theoretical analysis, which has recently been estimated to be 
AT sat , 2 = -370 ± 120 MeV dchen et alJl2009h . 


Next, we study the coupling constants for hyperons. The coupling constants for the_cr 
meson are chosen so as to g i ve reasonable hyper on potentials (IMivatsu et al.l l2013al. 12014c 


Schaffner fe Mishustinl Il996l: lYang fe Shenl 120081) . Using the baryon self-energy given in 


Equations (fT5l) and (fT6lh the potential for hyperon Y em bedded in the matter of baryon B 
is expressed by the Schrodinger-equivalent form (Jaminon et al. 19811): 


U { Y B) = Yf - E ° Y H + 


2 AT 


\Yt H - E 


0H1 ■ 


(63) 


Y 


where E^ (0)H is the direct term of the baryon self-energy for the scalar (the time component 
of the vector) part. Then, we can determine the coupling constants, g a y , following the 
values of potential depth around no suggested from the experimental data of hyper nuclci: 
U[ N) = -28 MeV, U^ N) = +30 MeV and Ui N) = -18 MeV Jschaffner et al.l 1 199^ . The 
scalar strange coupling constants, q^*v, are restricted by the relation T+r ~ ' ~ 2Ui A ' > ~ 


2 U. 


flSchaffner et al.lll994J: iSchaffner fc Mishustinl 1 19961: 1 Yang fc Shen 


f/| A ' ) ~ —5 MeV which has been implied by the Nagara event fjTakahash i et al . 2001) 


2008b We here take 


Furthermore, the coupling constants for the vector mesons to hyperons, g u y, 9<t>Y, and 
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g p Y , are given by the following SU(3) relations: 


9 lo a — S'ws — 


9<t>\ — 9(j> E 


1 + -\/3z fan 6( 

— tan 

1 + \/3z tan 


-<7wJV, 


9 wS 


■fi'wAf, g^s 


1 — \/2)Z tan 9 V 
1 + \/Zz tan 9 V 
\f?>z + tan 6 V 
1 + \/3 z tan 6 


9uiN, 


-QujNi 


9 p n — ~g P T, — .9pE, S'pA — 0, 


(64) 

(65) 

( 66 ) 


with the values of 9 V and z given below Equation f[62|) ( Mivatsu et al.ll2013a . 2014 : R.iiken et al. 
20101 ). Through the above relations, those couplings are determined once the values of 
(g^N, 9 p n) are given. For the tensor couplings, we may again use the ratios, f^y/guY, 
Ur/g^Y, and f P Y/g P Y, suggested by the ESC08 model. Furthermore, the coupling of f n y is 
also chosen to be the value in the ESC08 model. 

In the present calculations, we consider two parameter sets of the coupling constants 
for hyperons: one is the original set suggested by the ESC08 model (hereafter we call it 
ESC08Y), and the other is the values based on naive SU(3) symmetry (SU(3)Y). I 11 both 
cases, we use the coupling values, g a N, g^N, g<f>N and g p n, for SU(3) symmetry in Table El 
Then, in SU(3)Y, using Equations (K7H) - (j55j) and the values of (r/ wA r, g pN ) in Table El the 
coupling constants for hyperons are calculated. In cont rast, in ESC08Y, we adopt the values 
of (guN,g P N) presented in Table IV of Rii ken e t al. (120101 ) in the calculation of coupling 
constants for hyperons. The couplings related to the hyperons are listed in Table El 


5.2. Neutron star properties 


The properties of neutron stars ar e, in general, estimated by solving the the T olman- 
Oppenheimer-Volkoff (TOV) equation (|Tolmanlll939l: lOppenheimer fe Volkofflll939h . Since 
the radius of a neutron star is remarkably sensitive to the EoS at very low densities, we use 
the EoS for nonuniform matter below Wr = 0.068 f m~ 3 , wh ere nuclei are taken into account 
using the Thomas-Fermi calculation dMivatsu et ahl 2013bl ). 


In Figure [D we illustrate the particle fractions for hadronic matter in ESC08Y and 
SU(3)Y. Although all octet baryons are considered in the present calculation, only ap¬ 
pears at 0.675 fm~ 3 in ESC08Y, whereas, in SU(3)Y, the A firstly appears at 0.475 frn -3 , 
followed by 5” at 0.515 fm -3 . We also note that the Fock contribution suppresses the ap¬ 
pearance of hyperons as compared with the RH calculation including the a* and (p hi SU(3) 
flavor symmetry (see Figure 1 in Mi vatsu et ah (I2013al lb Due to the effect of the strange 
mesons, the threshold density of the 5“ production in ESC08Y is higher than that in our 
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previous calculation (see Figure 6 in Mivatsu et al. f 2013bh ). Furthermore, the densities at 
which a neutron star reaches the maximum mass are almost the same in both cases. 


The chemical potentials for the n, A, S - and e~ in hadronic matter are presented in 
Figure [2j When the f3 equilibrium conditions given in Equations f!5Tl)- (f56]) are satisfied, 
hyperons can be generated. At low densities, the chemical potentials for the A and 5“ 
behaves similarly in both cases, because, as explained in Section 15.11 the hyperon coupling 
constants are determined so as to reproduce the observed, potential depths at n 0 . However, 
their behaviors in ESC08Y show some difference from those in SU(3)Y above % ~ 0.4 
fm -3 . In particular, the chemical potential for the A in ESC08Y is slightly higher than 
that in SU(3)Y at middle and high densities. As shown in Figured] in the density region 
below 1.4 fm -3 , the A disappears in ESC08Y, while it appears at 0.475 fm -3 in SU(3)Y. 
This is mainly because the repulsive force due to the c o and (j) mesons strongly affects the 
chemical potential for the hyperons through the self-energy of the time component given 
in Equation (TT6|) . As already mentioned in Table [5] the coupling constant g u y in ESC08Y 
is larger than that in SU(3)Y, and thus the absolute value of the hyperon self-energy of 
the time component becomes larger in ESC08Y as the density increases, especially at high 
densities. It leads to the larger p a. Then, the hyperon creation in ESC08Y is suppressed 
compared with that in SU(3)Y. In addition, the E hyperon does not appear in both cases, 
because the E-hyperon potential in nuclear matter around no, U^\ is chosen to be repulsive. 
Furthermore, it is found that the chemical potential for the nucleon at high densities affects 
the onset of the hyperon productions, and the difference between the chemical potentials in 
ESC08Y and SU(3)Y becomes about 160 MeV around 1.0 fm -3 , around which a neutron-star 
mass reaches the maximum value. We note that, compared with the usual RH results, the 
Fock contribution also affects the chemical potential of neutron at high densities and hinder s 
the hyperon creation (IMivatsu et al.l 120121: iKatavama et al.l 120121: IWhittenburv et al.l 120121 ). 


The meson fields are presented in Figure [3] The (j) meson considerably contributes to 
the baryon interactions even at low densities because of the mixing effect in SU(3) flavor 
symmetry. In ESC08Y, the a* meson emerges above the density at which the hyperon 
is created, because we assume that g a *n = 0 in the present calculation. In contrast, the a* 


— ---~————j-- — ^ — v,- — 

reproduced even if the attractive force due to the a* is not included ( 

Schaffner et al. 

1994; 

Schaffner & Mishustin 

1996; 

Takahashi et al. 

2001 - 

Yang & Shen 

2008 

)• 


In Figures [4] and [5] we respectively show the EoS and the mass of a neutron star as 
a function of the neutron-star radius in ESC08Y and SU(3)Y. As well known, the inclu¬ 
sion of hyperons makes the EoS soft, and thus the maximum mass of a neutron star is 
reduced (see the curves for SU(3)Y in Figures [4] and [5]). However, due to the suppression 
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of the hyperon appearance and the large chemical potential for neutrons at high density, 
the EoS in ESC08Y keeps stiffness and the maximum mass of a neutron star can satisfy 
the recent 2constraint from the measurem ents of PSR J1614-2230 and PSR J0348+0432 
( Demorest et al. 201 q[ Antoniadis et al. 2013 1. Compared with our previous calculation de¬ 
noted by MYN13 (Miv atsu et al. 2013b!), the maximum mass in ESC08Y becomes heavier 
from 1.951M 0 to 2.029 M®. This is mainly because of the additional, repulsive force due 
to the strange, vector 0 meson. This fact can also be seen in th e result based on the RH 


calculation in SU(3) fla vor symmetry flMivatsu et al.l l2013al. I20l4 IWeissenborn et al.l 12011 


Lopes fe Menezesl l2014h In contrast, the maximum mass in SU(3)Y is below the 2M 0 con¬ 
straint from the pulsar observations. The radius of a neutron star around the maximum 
mass lies in the shaded areas obtained from theoretical analysis using the Markov chain 
Monte Carlo within Bayesian framework (Steiner et al. 2010), but the radius of the neutron 
star with 1.4M 0 in the present calculations is about 1 km larger than the range predicted 
by Bayesian analysis. 


5.3. Quark coexistence in neutron stars 

As shown in Figure (2] the chemical potentials for baryons increase as the density grows. 
Therefore, quarks may also be generated in the core region if the j3 equilibrium conditions 
given in Equations — (T5T]) are satisfied at some critical baryon density, In the MP 
where quarks coexist with hyperons as well as nucleons and leptons, we impose the conditions 
of equal pressure and global charge neutrality (see Equations (T58|) and (159]P in the range of 
the volume fraction, 0 < X < 1. 

In the present calculations, the phase transition from hadrons to quarks can be controlled 
by the density-dependent bag constant in Equation (1471) . In order to study the influence of 
the phase transition on the properties of neutron stars with hadrons, leptons, and quarks 
(i.e. hybrid stars), as shown in Table [6j we examine six cases, where the parameter j3 varies 
between 0 to 0.2, in ESC08Y and SU(3)Y. In Figure^)] the density-dependent bag constant is 
presented as a function of the total baryon density. It is found that the critical bag constant, 
B^ c \ and the critical chemical potential, //^, become smaller as the parameter 0 is larger. 

Relevant particle fractions for hybrid-star matter are given in Figures [7] and [HJ With the 
increase of the parameter 0 in the density-dependent bag constant, the threshold densities of 
quark productions move toward lower densities, and the hyperon population at high densities 
is suppressed in both hadronic cases, ESC08Y and SU(3)Y. Due to the quark productions, 
the S - completely disappears at the densities below 1.4 frn -3 , while, only in SU(3)Y, the 
A can still be generated even in the MP because of its neutral charge (see the panel (6) in 
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Figure [8]). In addition, the lepton fractions dwindle immediately at the densities above the 
onset of the quarks, because negatively charged particles are replaced by an abundance of 
quarks to conserve the charge neutrality. 

In Figure |9l we present the EoS for hybrid-star matter in ESC08Y. In general, the quark 
productions soften the EoS, as well known. Increasing the parameter /3, the EoS becomes 
softer because the threshold densities of the quarks move toward lower densities, as seen in 
Table 13 and Figure [0 Furthermore, the pure QP emerges at lower energy densities as the 
parameter f3 increases. 

We summarize the properties of neutron stars in Table [0 and illustrate the mass-radius 
relation of hybrid stars in Figure [10J It is found that, except the case (1), the maximum mass 
of a neutron star is reduced because the quark productions make the EoS soft. Moreover, 
the maximum mass is realized in the MP, and the transition to pure quark matter occurs 
only in the neutron star which already lies on the gravitationally unstable branch of the 
stellar sequence. Meanwhile, in the case (1), the threshold density of the quark productions, 
rig , follows the central density, n c , at the maximum-mass point, and then the maximum 
mass appears before the QP transition. In the case (2) of ESC08Y, the maximum mass 
reaches 2.OO3M 0 even if quark matter appears in the core of a neutron star. Therefore, 
the present EoS for hybrid-star matter can reasonably explain the recent mass constraint 
from astrophysical observations even if hyperons and quarks are included. As the properties 
of hybrid stars in SU(3)Y are similar to those in ESC08Y, we here do not present them. 
However, in SU(3)Y including quark degrees of freedom it is impossible to support the 
massive neutron stars (see Table ED- 


The neutron-star mass as a function of the total baryon density in ESC08Y is presented 
in Figure EU We find that, in spite of varying the parameter /3 in the density-dependent 
bag constant, all the maximum masses exist around the density of 1.0 fm -3 . In contrast, 
the threshold densities of the quark productions are extremely sensitive to the parameter /3, 
as already seen in Figures El and [HI and the end points of the MP also strongly depend on 
the parameter /3. Even if the quark degrees of freedom as well as hadrons and leptons are 
considered in the core of a neutron star, the EoS with the parameter range, /3 < 0.025, is 
consistent with the observation of heavy (~ 2M 0 ) neutron stars. From Table |H1 we can finally 
find the lower limit of the critical chemical potential for f3 = 0.025, /j, b ~ 1.5 GeV, which 
is almost the same value as in the calculation using the Polyakov-loop extended Nambu- 
Jona-Lasinio model with the vector-type four-b ody inte ract ion and the hadron resonance 
gas model with the volume-exclusion effect flSasaki et ah 2013 ). 
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6. Summary 


We have constructed the EoS for neutron stars with hyperons and quarks, and studied 
the properties of neutron stars with and without deconhned quarks in the core region. Under 
RHF approximation, the EoS for hadronic matter has been calculated by including of the 
strange (a* and </>) mesons as well as the light non-strange (a, u>, 7r, and p) mesons. In 
addition, we have used the chiral quark-meson cou pling (CQMC) model to take the variation 


of the in-m edium baryon structure into account (iNagai et al. 2008 : Mivatsu fe Saito 2009 


Saito 201 0). The EoS for quark matter has been calculated using the MIT bag model with 
one-gluon-exchange interaction, and the phase transition from hadronic matter to quark 


matter has been achieved under Gibbs criteria for chemical equilibrium (IGlendennind 1992 

200 lh . 


In the present calculations for hadronic EoS, we have examined the extension from 
SU( 6 ) spin-flavor symmetry based on the quark model to SU(3) flavor symmetry in deter- 
mini ng the isoscalar, vector-meson couplings to the octe t baryons within RHF a pproxima- 
tion ( Mivatsu et al. 2013a . 2014 : Weissenborn et al. 2011 : Lopes fe Menezes 2014). We have 
found that, in both cases of the SU( 6 ) and the SU(3) symmetries, the calculated proper¬ 
ties of nuclear matter at no, for instance, the incompressibility, the slope parameter of the 
symmetry energy, etc. are consistent with the experimental data and/or values predicted by 
other theoretical calculations. 

For hyperon coupling constants, we have adopted two parameter sets: one is the original 
set suggested by the ESC08 model (ESC08Y), and the other is based on naive SU(3) sym¬ 


metry (SU(3)Y). We have found that, compared with our previous results (IMivatsu et al. 


2013b), the maximum mass in ESC08Y becomes heavier from 1.95M 0 to 2.O3M 0 , because of 
the additional, repulsive force due to the strange vector ( 0 ) meson and the hyperon suppres¬ 
sion due to the large chemical potential of neutrons caused by the Fock contribution. While, 
in SU(3)Y, it is difficult to satisfy the 2M 0 constraint from the recent pulsar observations, 
because the hyperon appearance move toward lower densities, due to the small chemical 
potentials of the A and S _ , and it consequently makes the EoS softer. 


In order to study the effect of the phase transition on the proper 


have introduced the density-dependent bag constant for quark matter ([Burgio et. al, 


ies of hybrid stars, we 

2002 ®). 


If we presume the first order phase transition under Gibbs criteria ( Glendenning 19921" 2001 ). 


quark matter suppresses the population of hadrons and leptons, and it especially hinders the 
hyperon production. We have also found that the maximum mass in the case (2) of ESC08Y 
can reaches 2.003M© (see Table [TJ) even when hyperon and quark degrees of freedom are 
taken into account. We note that, even if quarks interact without strong, repulsive vector 
interaction, the present result can satisfy the mass constraint from the recent astrophysical 
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observations (jDemorest et aJj[2010|; lAntoniadis et al.l 120131) . Then, we have determined the 
lower bounds of the critical density for the phase transition, rtl? ~ 0.82 fm -3 , and the 


critical chemica l pot ential, // r ^ ~ 1.5 GeV, which are consistent with the recent calculation 
by Sasaki et ah ( 2013 ). In addition, our results have shown that the transition to pure quark 
matter occurs only in the neutron star lying on the gravitationally unstable branch of the 
stellar sequence. 

Finally, we give several comments on the future work. It may be important to consider 
the mixing effect of AS channel in the Fock diagram through the 7r or p meson exchange. 
Furthermore, the K or K* meson exchange between two baryons may be desirable in order 
to take the mixing of N-A, A-S, AT-E, or A-S into account, because these mixing may affect 
the EoS for hadronic matter, and con s equently the particle fr actions in the core region may 
be changed (jKatavama fe Saitol 120151: lYa/mamoto et al.l 120141) . For the possibility of other 
exotic degrees of freedom, it is als o interesting to include the effect of meson condensates 
flRvu et ahl 20071 ; Muto e t al.l 20091 ) in the RHF calculation for neutron stars. 


This work was supported by the National Research Foundation of Korea (Grant No. 
NRF-2014R1A2A2A05003548). 
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Fig. 1.— Particle fractions, Y, for hadronic matter in ESC08Y and SU(3)Y, which is defined 
as Yi = pi/riB with pi being the number densities of particle species i = B,£. The upper 
panel (a) is for the case of ESC08Y, and the lower panel (b) is for the case of SU(3)Y. The 
thick line shows the density at which a neutron star reaches the maximum-mass point by 
solving the TOY equation. 
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Fig. 2.— Chemical potentials for the n, A, H and e in ESC08Y and SU(3)Y. The filled 
circle denotes the onset of A or B“. The labels (a) and (b) are the same as in Figure [D 
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n B (fm' 3 ) 

Fig. 3.— Meson fields in ESC08Y and SU(3)Y. The labels (a) and (b) are the same as in 
Figure [D 
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Fig. 4.— Equations of state for hadronic matter in ESC08Y and SU(3)Y. We also present 
our previous EoS denoted by MYN13 (IMivatsu et al.1 2013bj). The Elled circle shows the 
point at which a neutron star reaches the maximum mass. 
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Fig. 5.— Mass-radius relations for hadronic stars in ESC08Y and SU(3)Y. The filled circle 
shows the maximum-mass point of a neutron star. The shaded bands show the recen t 2 
constraint from the measurem ents of PSR J1614-2230 and PSR J0348+0432 ( Demorest et al. 
2010l; lAntoniadis et al.l 120131) . The two shaded areas show the la and 2 a confidence limits 
for neutron-star radii of fixed mass obtained from the analysis by Steiner et ah (2 010 ). 
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Fig. 6.— Density-dependence of bag constant given in Equation fl47lh The circles (squares) 
show the critical points in ESC08Y (SU(3)Y) (see also Table E|). 
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Fig. 7.— Particle fractions, 17, for hybrid-star matter in ESC08Y. The labels (2)-(5) desig¬ 
nate the case number in Table |6j 



Fig. 8.— Particle fractions, 17, for hybrid-star matter in SU(3)Y. The labels (3)-(6) designate 
the case number in Table |6j 
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Fig. 9.— Equations of state for hybrid-star matter in ESC08Y. The filled circles show 
the points at the maximum mass of a neutron star, and the squares (triangles) mean the 
beginning (end) points of the mixed phase, i.e. y — 0 and y = 1 in Equations (16(1 and (j6Tj). 
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Fig. 10.— Mass-radius relations for hybrid stars in ESC08Y. The shaded bands and areas 
are the same as in Figure [5j The filled circles, squares, and triangles are explained in the 
caption of Figure EB 
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Fig. 11.— Mass-density relations for hybrid stars in ESC08Y. The shaded bands and areas 
are the same as in Figure [5j The filled circles, squares, and triangles are explained in the 
caption of Figure EB 















Table 1. Values of a B , b B , a' B , and b' B for the octet baryons in the CQMC model. 


B 

a B (fm) 

b B 

a' B (fm) 

// 

° B 

N 

0.118 

m 



A 

0.122 

Hi 

0.290 

1.00 

E 

0.184 

mm 

0.277 

1.15 


0.181 

1.15 

0.292 

1.04 


Note. — We assume that the scalar 
strange (a*) meson does not couple to 
the nucleon. 


Tabic 2. Functions Ai, Bi, Ci , and Di. 


i 

Ai 

Bi 

Ci 

Di 

a 

9l B {°)®* 

9l B (°)@* 

- 2 9l B (°)$* 


a* 

9l* B (v*)@a* 

S?-b(O 0 <7* 

~2 gl* B {v*)<&a* 


LOW 

2 9u B ®u 

^9u B ®u 

~~^9u B ^u 


LOTT 

- (fu B /2M) 2 ml<d u 

-3 (f uB /2M) 2 ml@ u 

4(W2 7W) 2 ^ 


loyt 




12 ( fu B 9u B /2M ) IV 

4>vv 

2 9 2 ^ B ®<j> 

~^9^) B ®(t> 

—^94> B ^0 


<f>TT 

- (Ub/ZM) 2 

-3(f<j> B /2M) 2 rri^O^ 

HUb/2 7W) 2 ^ 


4>vt 




12 {f<j, B g^ B / 2 M) r 0 

Pvv 

2 9p B ®p 

~^9p B @p 

~^9 p b^p 

... . 

Ptt 

-( f pB /2M) 2 m 2 p Q p 

-3(f pB /2M) 2 m 2 p Q p 

4(f pB /2M) 2 * p 


Pvt 




12 (/ pB g pB / 2 M) r p 

7tpv 

~fl B ^ 


2(f nB /m n ) 2 U 7T 



Note. — The index i is specified in the left column, where V(T) stands for the vector (tensor) 
coupling at each meson-BB' vertex. The last row is for the (pseudovector) pion contribution. 
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Table 3. Coupling constants related to the nucleon. 


Symmetry 

sIn/^ 

£jv/ 47r 

9%nI^ 

9 2 p nI^ 

SU(6) 

3.78 

5.74 


0.350 

SU(3) 

3.38 

5.41 

0.631 

0.412 


Note. — The (j) meson does not couple to the nucleon 
(. N ) in SU(6) symmetry, while it couples to N in SU(3) 
symmetry. 


Table 4. Properties of symmetric nuclear matter at tiq. 



m* n /m n 

I<o 

Jo 

-E'sym 

Lo 

K 

-* Y sym 

K 

-* Y asy 

-^■sat,2 

Symmetry 


(MeV) 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

SU(6) 

0.742 

275 

-366 

32.5 

75.3 

-48.0 

-500 

-400 

SU(3) 

0.747 

269 

-364 

32.5 

78.0 

-44.5 

-513 

-407 


Note. The effective mass of nucleon, incompressibility, third-order incompressibility, and 
symmetry energy are respectively denoted by Mjy, Kq, Jq, and E sym . The slope parameter and 
curvature parameter of the symmetry energy, Lq and A" sym , are also listed. Using the parabolic 
approximation for the EoS, the 2nd derivative of the isobaric incompressibilit y co efficient is 
fru7Dn r ^ — A T ^ find nQTQinott 


Aerawal et al. 

A _ 

2012; 

l psy- Ko^u “ - " 

Sulaksono & Aerawal 

2012) 


asy 


l sym 
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Table 5. Coupling constants related to the hyperons. 


Set 

ESC08Y 

SU(3)Y 

9a A 

2.48 

1.74 

9aT, 

1.98 

1.17 

[I] 

b 

1.95 

1.35 

9a* A 

0.876 

# a 

9a*?, 

0.876 

. a 

[I] 

b 

1.05 

# a 

9uiA 

2.82 

1.85 

9u? 

2.82 

1.85 

£ 

[I] 

2.09 

1.37 

9<$>A 

-2.16 

-1.42 

9<t>? 

-2.16 

-1.42 

9<t> H 

-3.11 

-2.04 

9p? 

1.38 

1.28 

9pS 

0.692 

0.642 


Note. — All the values are divided by \[Am. We list the original set 
suggested by the ESC08 model (ESC08Y) and the calculated set based 
on naive SU(3) symmetry (SU(3)Y). We assume g a *\ = and g p a = 0. 

a Because the a-meson contribution in SU(3)Y already gives U^ ) = —19 
MeV and U ( £ } = —15 MeV at n o, the additional, attractive force due to 
the a* meson is not required. 







Table 6. Phase transition properties at the critical density, (fm 3 ). 


Case (f3) 

ESC08Y 

SU(3)Y 

( c ) 

n B 

£(c) 

( c ) 

( c ) 

n B 

£(c) 

( c ) 

/4 

(1) 0.000 

1.065 

400 

1702 

1.775 

400 

2031 

(2) 0.025 

0.820 

224 

1474 

1.185 

131 

1605 

(3) 0.050 

0.680 

184 

1350 

0.945 

105 

1444 

(4) 0.100 

0.565 

143 

1246 

0.660 

107 

1274 

(5) 0.150 

0.500 

123 

1188 

0.505 

121 

1188 

(6) 0.200 

0.450 

115 

1147 

0.450 

115 

1147 


Note. — We list six cases, where the parameter j3 varies between 0 to 0.2, for 
ESC08Y and SU(3)Y. The critical bag constant, B^ c \ and the critical chemical 
potential, , which are in the unit of MeV, are calculated by Equation (1471) . 


Table 7. Properties of a neutron star in ESC08Y and SU(3)Y. 


Case (/ 3) 


ESC08Y 



SU(3)Y 


■^max 

M ma , x / Mq 

n c 

■^max 

M max /M 0 

n c 

(1) 0.000 

11.13 

2.029 

1.035 

11.36 

1.874 

1.020 

(2) 0.025 

11.56 

2.003 

1.005 

11.36 

1.874 

1.020 

(3) 0.050 

11.71 

1.958 

0.975 

11.46 

1.873 

1.025 

(4) 0.100 

11.63 

1.896 

1.010 

11.51 

1.863 

1.030 

(5) 0.150 

11.50 

1.866 

1.050 

11.43 

1.856 

1.060 

(6) 0.200 

11.41 

1.853 

1.070 

11.37 

1.850 

1.075 


Note. — We list the neutron-star radius, -R max (in km), the ratio of the neutron- 
star mass to the solar mass, M max /M 0 , and the central density, n c (in fm -3 ), at 
the maximum-mass point. 














